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Edge states of semi-infinite nanowires in tight binding limit are examined. We argue that under- 
standing these edge states provides a pathway to generic comprehension of surface states in many 
semi-infinite physical systems. It is shown that the edge states occur within the gaps of the corre- 
sponding bulk spectrum (thus also called the midgap states). More importantly, we show that the 
presence of these midgap states reflects an underlying generalized supersymmetry. This supersym- 
metric structure is a generalized rotational symmetry among sublattices and results in a universal 
tendency: all midgap states tend to vanish with periods commensurate with the underlying lattice. 
Based on our formulation, we propose a structure with superlattice in hopping to control the number 
of localized electronic states occurring at the ends of the nanowires. Other implications are also 
discussed. In particular, it is shown that the ordinarily recognized impurity states can be viewed as 
disguised midgap states. 

PACS numbers: 73.20.-r, 73.20.At, 73.21.Hb 



I. INTRODUCTION 

The one-dimensional (ID) wire has been of great the- 
oretical and experimental interests in the past. This is 
because of not only the wide variety of fascinating phe- 
nomena it exhibits, but also the testing grounds it offers 
for ideas that may become applicable in higher dimen- 
sions. In practice, ID wires needs not be physically one 
dimension. It may result from projection after a partial 
Fourier transformation from higher dimensional models. 
For instance, a superlattice structure can be reduced to 
an equivalent ID structure after a partial Fourier trans- 
formation along the direction normal to the layers. Sim- 
ilar examples include d-wave superconductors, graphite 
sheet, and many other systems. Therefore, understand- 
ing the ID wire is an ideal first step toward the un- 
derstanding of any higher dimensional problems. Fur- 
ther boost for studying ID wires comes from recent ad- 
vances achieved in nanotechnology. Here the feasibility 
for bottom-up assembly of single nanowiresi has made 
direct investigation of finite ID wires possible. Never- 
theless, conventional studies of the ID wire have mostly 
been focused on its bulk properties, whereas assembled 
nanowires can only have finite lengths and must termi- 
nate at some sites (the ends, or the edges). It is therefore 
desirable to reconsider the effects of the ends to the prop- 
erties of the nanowires. 

The commonly recognized edge (or surface) effects in 
the physics of nanostructures are concerned with the 
large volume fraction of the boundaries. However, from 
either fundamental or practical viewpoints, the possible 
occurrence of edge modes and its influences on the prop- 
erties of the system poses a much more interesting prob- 
lem. For example, when applying carbon nanotubes as 
emitters for screen displays, the occurrence of edge states 
may change the density of electrons at the edge and thus 
affects the threshold working potential. It is therefore 



of great technological interest if one could devise a way 
to engineer the number of edge states. From the funda- 
mental viewpoint, the elegant role of edge excitations in 
the physics of Quantum Hall systems^ is a well known 
example that illustrates the importance of edge states. 
Generally speaking, the edge states occur within the gap 
of the bulk energy spectrum and are called the midgap 
states. The existence of these states causes anomalous 
properties near the end which can manifest in tunnel- 
ing measurements. A recently discovered example is the 
zero-bias conductance peak observed in the dl /dV mea- 
surement of the metal-G?-wave superconductor junctions 3 . 
When electron-electron interactions are present, as oc- 
curs for an externally-implemented magnetic impurity, 
"intrinsic" Kondo effects may also arise due to these lo- 
calized states, causing zero-bias anomaly near the Fermi 
energies—. Furthermore, if the system is finite, coupling 
between edge states can not be neglected. An example 
is the anomalous paramagnetic behavior observed in car- 
bon nanoribbons^, where we have recently shown that 
there are residual antiferromagnetic couplings between 
edge spins in this system^. All these examples clearly 
illustrate the important role of edge states in the physics 
and applications of nanostructures. 

In previous work, applying the Green's function ap- 
proach, we have shown that broken reflection symme- 
try is a necessary condition for the occurrence of edge 
states, and the energies of edge states are the roots to 
the Green's function 3 .. In this work, resorting to the su- 
persymmetric method, we further develop a systematic 
way to determine the wavefunctions and the precise en- 
ergies of the edge states. 

Conventionally, the usage of the supersymmetric 
method in the condensed matter physics has been fo- 
cused on applying the supersymmetry (SUSY) quantum 
field theory to disorder systems^. The application of the 
corresponding (0+1) dimensional limit -the SUSY quan- 
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turn mechanics 8 , however, is quite limited. Nevertheless, 
it has been realized^ that the zero-bias anomaly in d- 
wave superconductors is closely related with the SUSY 
quantum mechanics. These studies, however, are done in 
in the continuum limit, using the semi-classical approxi- 
mation, while the more relevant limit for high Tc super- 
conductors and many other systems is the tight binding 
limit. Furthermore, the zero-energy state was the pri- 
mary focus, while not all the states localized at the edge 
have zero energy. It is therefore important to see if the 
idea of SUSY quantum mechanics can be generalized to 
understand the finite-energy midgap states, in particular 
those in the discrete condensed matter systems, as well. 
In this work, we shall show that indeed, this is possible. 
We shall first show that the semi-infinite tight-binding 
d-wave superconductors belongs to a more general class, 
the bipartite system, and which can be well described by 
the conventional SUSY quantum mechanics^. Here the 
supersymmetric partners are two sublattices of the same 
system and the SUSY is characterized by a hermitian su- 
percharge Q and the SUSY Hamiltonian H s = {Q, Q}/2 
with [H S ,Q] = 0. For bipartite systems with nearest 
neighbor hoppings, Q is identical to the physical Hamil- 
tonian (= H2) and hence H s is a quadratic functional 
of Hi . Furthermore, the zero-energy state is annihilated 
by the supercharge, which then constitutes one of the 
conditions for determining the zero-energy state; while 
the other condition is to require it to decay from the 
edge. It is found that this conventional SUSY quan- 
tum mechanics can be appropriately extended to describe 
the semi-infinite p-partite systems with nearest neighbor 
hoppings. First, when p > 3, the original supercharge 
splits into two: In addition to the physical Hamiltonian 
Hp, a second supercharge Q p can be formed. They both 
commute with the SUSY Hamiltonian Hs- Furthermore, 
only when p — 2, Hs ~ H s is a quadratic functional 
of Hp. In general, Hs is a polynomial functional of 
Hp. This is a reminisce of the fractional SUSY quantum 
mechanics^ in which the SUSY Hamiltonian is general- 
ized to be integer power of the supercharge. Neverthe- 
less, our model is different and provides more realistic 
generalization of the conventional SUSY. The upshot of 
this generalization shows that, in addition to the zero- 
energy state, all the midgap states, including finite en- 
ergy ones, are annihilated by the supercharge Q p . The 
wavefunctions of the midgap states thus obtained tend 
to vanish with the same period commensurate with p: 
vp ps (• • •, 0, • • •, 0, • • •, 0, • • •). These zeroes cut the orig- 
inal Hamiltonian into smaller ones so that the energies 
of the midgap states are determined by the eigenvalues 
of the Hamiltonian within each period. As a result, the 
matrix for determining the energies of midgap states is 
of size much smaller than the size of the original Hamil- 
tonian. This reduction in matrix size heavily reduces 
the computation for determining the occurrence of the 
midgap states and provides a way to control the occur- 
rence of the midgap states. As an application, we propose 
a structure with superlattice in hopping with period p to 



control the number of localized electronic states occur- 
ring at the end of nanowires. In that case, the number 
of edge states is simply p — 1. 

As the period p goes to infinity, the ensemble of config- 
urations of hopping forms a semi-infinite disorder chain. 
This limit has been extensively investigated during the 
pasliii since Dyson's seminal workiS in which it was 
pointed out that the average density of state (DOS) is 
enhanced at zero energy. From our point of view, this 
enhancement also reflects that the system has high prob- 
abilities to take the above-mentioned form for the ground 
state. The presence of the boundary breaks translational 
invariance. Thus, unlike the bulk case where the DOS at 
zero energy has no spatial dependence, the enhance DOS 
at zero energy for semi-infinite disordered wires has the 
largest amplitude near the edge. Even for slight disor- 
ders, the effects of enhanced DOS at zero-energy are still 
observable. This offers a possible explanation for many 
unexpected zero-bias anomalies observed in tunneling ex- 
periments because, unless extremely carefully controlled, 
junction qualities are usually rather poor and disorders 
can easily set in near the junctions^. 

Other implications and extensions of our generalized 
SUSY quantum mechanics will also be discussed. In par- 
ticular, we shall demonstrate that by appropriate map- 
pings, the ordinarily recognized impurity state can be 
viewed as a disguised midgap state. Such mapping pro- 
vides a simple way to construct the impurity wavefunc- 
tion and the corresponding energy. In addition to this 
application, possible extension to include the electron- 
electron interactions will also be discussed at the end of 
this paper. 

This paper is organized as follows. In Sec|H] we lay 
down the basic tight binding model considered in this 
work and illustrate the SUSY quantum mechanics for 
the bipartite systems. In Sec llIII we generalize the super- 
charge and supersymmetric Hamiltonian to the p-partite 
systems and discuss the disorder limit. We also point it 
out of how to engineer the number of edge states by using 
a superlattice structure. By applying the SUSY quantum 
mechanics, we illustrate in Sec II VI how an impurity state 
can be viewed as a midgap state. In Sec[V] we conclude 
and discuss possible generalization to include electron- 
electron interactions. Appendices A and B are devoted 
to technical details of superalgebra and computation of 
commutators. 



II. THEORETICAL FORMULATION AND 
SUPERSYMMETRIC QUANTUM MECHANICS 

We start by considering the ID atomic chain as illus- 
trated in Fig. da). This is the most general ID atomic 
chain in which reflection symmetry with respect to the 
edge point is broken and, consequently, edge states might 
arisA In the tight-binding limit, we consider the follow- 
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ing Hamiltonian to model this system: 
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Here the subscript p indicates the period of the lattice 
and i is the site index; tj is the hopping amplitude be- 
tween site i and its nearest neighbors, c; (cj) is the elec- 
tron annihilation (creation) operator, and Vi is the lo- 
cal potential at site i. We shall assume that both 
and Vi are periodic with period p, namely t p+ i = ti 
and Vp+i = Vi. In real systems, this Hamiltonian may 
correspond to an assembly of p different atoms repeat- 
edly arranged into a line (see Fig. For wires com- 
posed of atoms of a single species, H p may describe sys- 
tems which exhibit density- wave order. This includes 
polyacetylene 14 , which has a dimerized structure and 
corresponds to p — 2, and polymers with higher com- 
mensurability charge density wavesit In the following, 
we shall call p = 2 the t\-t2 model, and similarly for 
models with higher periods. As mentioned in the intro- 
duction, Hp may also represent the reduced model of a 
higher dimensional structure after partial Fourier trans- 
formation. For example, for a semi-infinite graphite sheet 
with zig-zag edge, since the system is translationally in- 
variant along the edge, a partial Fourier transformation 
can be applied along this direction, leading to an effec- 
tive ID model. In this case, it is identical to the t\- 
t% model except that now t\ and t-i are fc-dependcnt 3 : 
t\ = 2io cos (V3 fey a/2), £ 2 = to , where a is the lattice con- 
stant and k y represents the Fourier mode. This approach 
has been successfully applied to understand the anoma- 
lous properties near the edge in carbon ribbons^. As a 
final example, we note that the operator c$ in H p needs 
not be restricted to be the electron annihilation operator. 
For example, after applying the Jordan- Wigner transfor- 
mation, one can map a ID quantum XY spin chain to 
a ID model described by H p . Specifically, we have U 
replaced by the exchange coupling for nearest neighbors 
Ji/2, and vi replaced by the local magnetic field hi. It is 
clear from these examples that H p is quite general and 
captures the physics of many interesting systems. 

To investigate the behavior of H p near the edge, as 
a first step, we calculate the local density of states at 
the end point using the generalized method of image de- 
veloped in Ref. 3. Fig. |2] shows typical local density of 
states at the end point for small periods. The param- 
eters are carefully chosen so that all possible midgap 
states are present. In particular, we have set Vi = 0, 
which amounts to choosing the energy zero as the ori- 
gin. These results show that midgap states are indeed 
the most prominent features at the end point. To under- 
stand how the midgap states arise, we first investigate 
the ti-t% model with m — in details. In this case, since 
the lattice is bipartite, it is convenient to distinguish the 
amplitudes at the odd and the even sites by writing the 
wavefunction as \& = (0 O , </> e ). The Hamiltonian then 
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FIG. 1: Schematic plot of (a) an assembled atomic chain and 
(b) the corresponding models with small periods: p = 2, 3, 
and 4. 




FIG. 2: The local density of states at the end point for small 
periods: p = 2, 3, and 4. The parameters are carefully chosen 
so that the midgap states are manifested. Top: ti = 2.0, 
ts = 3.0 = 2); Middle: ti = 2.5, t 2 = 1.8, t 3 = 3.0 (p = 3); 
Bottom: ti = 2.4, ta = 2.8, t 3 = 1.3, U = 3.4 (p = 4). All 
potential Vi = 0, and we have included a lifetime 8 = 0.02. 
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Here is the null matrix and A is a non-Hcrmitian ma- 
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It is interesting to note that the adjoint of A satisfies 
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Here the operator e effectively reflects the wavefunction 
with respect to the mid point of the lattice. 

In the case of infinite chains, it is not hard to check 
that the corresponding matrices A and A commute with 
each other and hence can be diagonalized simultaneously 
in Fourier space. For semi-infinite chains, however, A 
and A do not commute and the spectrum of the t\-t2 
model can be best understood in terms of the supersym- 
metric quantum mechanics 8 . For this purpose, we first 
identify H2 as the supercharge Q2, which connects even 
and odd sites. The block-diagonal matrix (H2) 2 (= H s ) 
is then identified as (up to a factor of two) the corre- 
sponding supersymmetric Hamiltonian, whose diagonal 

blocks = AA and ijf = A* A are, respectively, 
the effective Hamiltonians for the odd and the even sites. 
Note that because A and A do not commute, H^ ^ fff . 
We will show below that the difference between Hg and 
-fff is the origin of the midgap states. Obviously H s is 
positive definite with the possibility when its spectrum 
touches zero. When the later happens, the ground state 
energy of H s vanishes, the ground state wavefunctions 
4> e and 4> [*o — (<^oj <Pe) ] have to be the zero-energy 
eigenfunction of A and A , ie. A0 e = and A <f> = 0. 
In other words, the supercharge annihilates the ground 
state wavefunction "to 

Q 2 *o= #2*0 = 0. (5) 

Clearly, in this case, the system has good supersymmetry 
because the ground state is invariant under "rotation" 
between even and odd sites: 

e««fetf = ¥ , (6) 

where 9 is any real number. The non-Hcrmiticity of A 
and A implies that forward and backward hopping am- 
plitudes between two sites are different, and hence the 
eigenf unctions have to either grow or decay from the end 
point. Obviously, because of the relation A + = eAe, any 
non-trivial eigenfunctions satisfy <j) e — £<p . Therefore, 
if 4> decays from the edge, <fr e must grow from the edge 
(vice versa). For semi-infinite chains, only the even sites 
are connected with the hard-wall boundary point. Thus 
4> e is forced to vanish while <f) decays into the bulk, so 
that ^0 = (0o,O). Note that the other possible state 
^0 = (0, 4> e ) resides on the other end of the chain and is 
pushed to infinity. Therefore, overall speaking, there is 
only half a chance for the existence of the ground state 
((/> o ,0). This also reflects in the hopping strength differ- 
ence. Indeed, we find that (j) decays only when t\ < <2- 
In this case, has a non-trivial zero energy eigenfunc- 
tion, while does not. Therefore, the system has good 
supersymmetry with the ground state ^0 being a local- 
ized state. For finite energies, however, <p e and <p need 
not be eigenfunctions of A and A . Nevertheless, the 
supersymmetry allows a simple and elegant way to find 
the whole spectrum for the case p = 2. This is because 
has the exact form as Hi (p = 1) with t{ — t\t2 
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FIG. 3: The effect of the potential is to break the particle-hole 
symmetry so that two side bands are distorted. However, the 
midgap state is not changed if v — 0. The parameters are: 
ti = 0.7 and ti = 2.0. For solid line, v = 0, v e = 0, while for 
dash line v = 0, v e = 0.5. 



and Vi = (t 2 + Since this is just the ordinary uniform 
hopping model, one can easily write down the eigenstate: 
4>e{n) = sin2nfc. The wavefunction at odd site can be 
then found by using the supercharge operator. We find 
that <fi — A(f> e /E with E being the spectrum of H% which 
satisfies £ 2 = t\ +t! + 2tit 2 cos2fc. Since E 2 > (ii~t 2 ) 2 , 
an energy gap opens up around E = when t\ ^ t%. In 
the case of t\ < ta, the ground state ^0 then arises as a 
midgap state. Note that H^ is almost identical to Hf 
except for the potential energy v\ = t\ at the end point; 
the deficit energy t\ is entirely due to the missing bond 
cut off by the boundary. We will elaborate on this in 
SecED 

We now address the effects of the potential Wj. For 
p = 2, it is convenient to denote the potentials over the 
even sites v e and the odd sites v . This decomposition, 
however, renders the particle-hole symmetry invalid at 
the level of the supercharge H 2 - Nonetheless, the spec- 
trum (E) of H2 can be mapped to the original spectrum 
of H s with Vi = (= Eg). For E s ^ this mapping is 
given by E s — (E — v e )(E — v a ), while for E s — 0, since 
<j) e = still holds, one has E — v a . Hence even though 
the physical spectrum E may have no particle-hole sym- 
metry, after appropriate transformations, the symmet- 
ric structure can be restored in E s ; in particular, the 
midgap state survives as clearly demonstrated in Fig. [3] 
For higher periods, the same manipulations as above can 
lead to similar conclusions. Therefore, unless explicitly 
needed, we shall ignore Vi in the following. 

Let us now apply the supersymmetric method to semi- 
infinite superconductors^. After partial Fourier trans- 
formation along the interface, the problem becomes ID 
superconductors with an end point. In this case, it is con- 
venient to write the wavefunction by "J = (u, v) with u = 
(til, U2, U3, • • • ) being particle- like and v = (vi, t>2,i>3, • • • ) 
being hole-like wavefunctions. The reduced mean-field 
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(BCS) Hamiltonian is Dirac-lika^ and can be generally 
written as 



Hbcs = 



M Q 

Q -M 



(7) 



where M corresponds to the reduced ID Hamiltonian for 
particles and Q is essentially the pairing potential. One 
can also rewrite Hbsc — M <£> a z + Q <g> a x , and treat 
this problem as a spin in the "magnetic field" (Q,0,M) 
pointing in the x — z plane. This analogy suggests that it 
is possible to rotate the magnetic field to the x — y plane. 
Indeed, this can be achieved by a rotation of 27r/3 with 
respect to the axis (1,1,1). The transformation matrix- 
is 



U 



1 

vi 



i i 

i — i 



(8) 



where 1 and i are semi-infinite matrices. The rotated 
Hamiltonian then takes the form of a supercharge like 
H 2 



H'r 



BCS 



BCS 



u = 



° M iQ ]. (9) 



M+iQ 



The wavefunction is rotated accordingly: <!>' = W^. 
Therefore, in the supersymmetric form, particles and 
holes are mixed. As an illustration, we consider the 
mean-field Hamiltonian for d-wave superconductors 

H R = toc^c^ + A^^Cji-c^Cj^+h.cilO) 

where (ij) denotes the nearest-neighbor bonds, to and 
Aij are, respectively, the corresponding hopping and d- 
wave pairing amplitudes. For the (1, 1,0) interface, after 
Fourier transformation along the interface (which is taken 
to be the y direction), we obtain^ 



A = M iQ 



-[i -t + d 
—t — d —ji —t + d 
-t-d -/i 



(11) 



where \i is the chemical potential, t = 2t cos(k y a/\/2) 
and d = 2Ao sm(k y a/^/2). The model for tight-binding 
d-wave superconductors is unique in the sense that the 
non-Hcrmiticity of A and A can be removed by a gauge 
transformation. For this purpose, we write t — d = te 9 
and t + d = te~ 9 with t = Vt 2 - d 2 and e 2s = (t - 
d)/(t + d). The eigenfuction of A is thus the gauge 
transformation of the eigenfuction 4>e of 
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-t -/j, 
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Specifically, we obtain v E'_E(n) 



-ng. 



(for A 



obtains ^ #(n) 



is). Furthermore, and <pE have 



the same eigenvalue E. This implies that the calculation 
of the zero mode is related to the spectrum of Ao. If 
we restrict our discussion to the propagating modes, the 
spectrum of Ao is simply the ordinary cosine band. We 
find that when — 2V t 2 — d 2 < /i < 2\Jt 2 — d 2 is satisfied, 
A and A can, respectively, support zero-modes of the 
form (0, e~ n9 (f>E=o) and {e n9 4iE=o, 0). In this case, the 
ground state of H s is selected by the sign of g, and thus 
the zero-energy midgap state is given by 



u(n) = —= \ j— 7 I sm(fcf)io), 



V2\t+\d\ 
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(13) 
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Here fei? is determined by —2y/t 2 — d 2 cos(fcj7fl) = A 4 and 
depends on k y ; therefore, even though the midgap states 
for different fcj/s have the same zero energy, their wave- 
functions have k y dependence. Note that for demonstra- 
tion, we have only considered the case when fcp is real. 
Complete solutions, however, require to include the situa- 
tion when fcp is complesi^. In both cases, the supersym- 
metry structure enables one to write down the explicit 
form of the zero-energy mode near the interface (1, 1, 0). 



III. GENERALIZED SUPERCHARGE AND ITS 
CONSEQUENCES 

We now generalize the above results to higher periods 
p > 3. It is useful to decompose the wavefunction as W = 
(4>i, 4>2,' • • ,4>p)i where <fi n denotes the sub- wavefunction 
formed by {^(kp + n); k = 0, 1, 2, . . . }. The Hamiltonian 
is then cast in the form 
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Again, here and A nm are block matrices; for all n ^ p, 
A nm = t n l are diagonal, while for (m,n) = (l,p) 
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To understand what happens for the semi-infinite chain, 
it is useful to start from the infinite chain with Hamil- 
tonian . In this case, also takes the same form 
except that A nm are further extended to i = — oo. If we 

remove the hopping strength t n and combine the remain- 

t 

ing A nm with A nm into Q nm for all m and n pairs, Q nm 
form a superalgebra if modulo p is performed (see Ap- 
pendix A for mathematical details). The energy bands 
of H£° are determined by 
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P(E, k) = det 
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where E is the energy and k is the Fourier mode. In gen- 
eral, there have at most p energy bands. However, since 
in the polynomial P(E, k), (— I) p+1 2t±t2t3 ■ ■ -t p cos(pk) is 
the only term that depends on k, the function P(E, k) — 
(— l) p+1 2tii 2 *3 • • • t p cos(pk) maps p bands into one sin- 
gle band: 2tit2t3 ■ ■ ■ t p cos(pk). This important ob- 
servation implies that when H p — the operator 
H s ee P(H p ,0) - (-l)P+ 1 2t 1 t 2 i 3 • • -t p is block-diagonal 
[there are p blocks with one for each fa, see Eas. l|B 1J) 
and l|B 2| l] and folds the spectrum of H^ into one sin- 
gle bandit. Therefore, Hs is similar to the supersym- 
metry Hamiltonian Hs- Indeed, for p — 2, we find 
H s = iff — {t\ + tf) which is essentially H s . 

For infinite chains, Hs is highly symmetric. In fact, it 
commutes with all Q m n- This reflects that it is symmet- 
ric under the permutation of n and m but it is more than 
that because any linear combination of ^2 tmnQmn also 
commutes with Hs- For semi-infinite chains, however, 
the above symmetry is broken: Not all Q mn commute 
with Hs- Physically, this is obvious because now fa is 
special and is the only component that connects with 
the boundary point i — directly. As a result, Hs is 
not completely block-diagonalized. In fact, because even 
for the infinite chains, fa(n) = sin(fcrtpa) is a solution for 
the pth block and it satisfies the hard-wall boundary con- 
dition at i = 0, the pth block is not affected. Therefore, 
there are only two blocks: one for the space formed by fa; 



the other mixes 



..and 



' P -i- 



This is demonstrated 



in Eq. l|B 1|> where we denote the block Hamiltonians by 
H+ and H- . 

Clearly, the t\-t2 model is special because H + and H 
are of the same size so that Hs is completely block- 
diagonal. This is where the usual SUSY quantum me- 
chanics applies. For p > 3, H + and H are not of the 
same size, a generalization of SUSY quantum mechanics 
is needed. First, it is important to see if one can find an 
operator, similar to the supercharge Q2, that commutes 
with Hs- H p is obviously a solution because H 2 = Q2- 
However, in analogy to the case of p = 2, a second super- 
charge by collecting all block matrices in H p that con- 
nects fa to other components can be formed: 












\A lp 



A-p—l,p 





for p > 3. (18) 



Note that the above defintion can also include p = 2. In 
that case, one squeezes the block Ai p into A12 to obtain 
the form of H2 in Eq.©. Because A12 = til and Ai p is 
given by Ea. Hl(j|) with t p being replaced by t 2l adding A12 
and Aip reproduces A defined in Eq.© precisely. Hence 
Ea. (|18|) can be regarded as an "analytical continuation" 
of Q2 to p > 3. Furthermore, as shown in Appendix 
B, [Hs,Q p ] = is satisfied, thus Q p provides a faithful 
generalization of Q 2 - It coincides with H p only when 
p = 2. Note that from Appendix B, one can actually 
see that out of Q mn contained in H p , Q p and H p (and 
their linear combinations) are the only two generators 
that commute with Hs for general p. 

Applying the condition that the supercharge annihi- 
lates midgap states <3/ci) one finds that 



and 



t 

^p—itprp- 



0. 



(19) 



<p p = implies that the wavefunction has the form 5*0 ~ 
(• • •, 0, • • •, 0, • • •, 0, • • •) as pointed out earlier—, while the 
second condition relates 4> p —i to <p±. It is important to 
realize that because Q p no longer coincides with H p for 
p > 3 , Qp alone does not determine the energies and 
wavefunctions of the midgap states. Instead, because 
Qp^o = 0, the operator H p — Q p determines the energies 



and further provides relations between fa, 



and 



fa . This analysis shows that the energies E m of midgap 
states can be different from zero and must satisfy 



det 



-E„ 

h 




V 



h 

-En 
t 2 







t-2 

-E n 



tp-2 



(20) 



tp-2 



-E m J 



Therefore, there are at most p — 1 midgap states. To 
stabilize the midgap states, one further requires *o to 
decay away from the edge. In the case of p = 2, this 
results in the condition t\ < t2- For p = 3, one first 
obtains from Ea. H2UI) E m = ±ti and fa — ±fa, which, 
when combined with Ea. (|19fl . results in ^(3k + n) — 
±t3/t2^{3k — 3 + n). Thus midgap states exit only when 
t2 < ta. In general, one needs to relate ipp-i to ip\. This 
further reduces the matrix in Eq. I|2(J|I and by defining the 
{p — 2) x [p — 2) matrix 



E m 


h 





h 


—E m 


t2 





t 2 


-E, 











(21) 



-E m } 
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we find that ipi = — p _ 2 V ! p-i' When combined 

with Ea. ll9fl . we obtain ipi — t p _2tp/t p -ihi 2 tpp+i. 
Hence the midgap state with energy E m exists only 
when i p _2^ph^p_ 2 > i p -i 19 . Note that for higher peri- 
ods, commensurate structures may appear in sublattices. 
These structures resemble the SUSY structures in lower 
periods. For example, when p = 4, there are at most 
three midgap states at E m — 0, ±\Ai +t\. In this case, 
the Hamiltonian H\ is already block diagonal in even 
and odd sites. For even sites, H\ is period of two and 
belongs to H 2 with t[ — M4, t' 2 = t%ts, fi[ = t\ + t 2 , and 
p! 2 = i| + t\ with the midgap energy given by E' m — fj,[ = 
t\ + t 2 . It is clear that E' m is precisely the square of 
E m = iy'if + if, and thus ±-\Ai + t\ are the midgap 
states resulting from the supersymmetry within the even 
sublattice. The remaining midgap state at E m = orig- 
inates from the same supersymmetric structure as that 
for p = 2 between even and odd sites. Hence <fi e vanishes 
for E m = 0. Obviously, the same effects may occur for 
any p that is not a prime number. 

As an application, we point out that a simple way to 
engineer localized edge states in a nanowire (such as car- 
bon nanotubes) is to introduce impurity-bonds periodi- 
cally and form a semi-infinite superlattice structure. In 
this case, one has t\ = t% = ■ ■ ■ = = t while t p = if . 
It is known that in the infinite case, there are p energy 
bands. However, when it becomes semi-infinite, p — 1 
localized edge states may arise simultaneously and they 
are located exactly in the p — 1 gaps among these en- 
ergy bands. The proposed superlattice structure thus 
provides a convenient way to control the number of edge 
states. The above result can be quite easily derived in 
our formulation. First of all, Ea. (|20|l implies that the 
energies of edge modes are exactly the energies of a finite 
atomic chain with uniform hopping t. In other words, 
E = 2t cos(fca) with ka = mn/p, m — 1,2, • • • ,p — 1. 
Note that these fc's occur exactly at the zone boundaries 
of the energy bands, hence they appear within the energy 
gaps. Ea. H19[) further implies that for any edge mode, 
their wavefunction satisfies ^(p + 1) = —t/t'^>(p — 1). 
Since ^(1) and — 1) have the same amplitude, we get 
|*(p+ 1)| = \t/t'\\^{l)\. Therefore, as long as \t/t'\ < 1, 
all midgap wavefunctions decay and hence the p — 1 
midgap states appear at the same time. 

Another example is to consider the limit when p goes 
to infinity. Since for any given configuration of {ti,i — 
1, 2, 3 . . . }, ti (1 < i < p = 00) is allowed to be any 
number; effectively, the wire is a disordered semi-infinite 
wires. It is well known that the density of states at zero- 
energy become enhancediii^ in the disordered wires. In 
the following, we will find that such enhancement can be 
also understood from the above point of view. Essen- 
tially, this is because for any given ti configuration, the 
system has high probability to settle into the form of the 
ground state wavefunction discussed above. First, be- 
cause the boundary breaks the translational invariance, 
if one decomposes the wavefunction as ^ = (<f> , (f> e ), it 
is still true that only the even sites are connected to the 
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-3-2-10 1 2 3 

Enerav 

FIG. 4: Averaged local density of states over 10000 samples 
at the end point for semi-infinite wires with disorders near the 
edge point. Only 5 lattice points are imposed with disorders. 
Here random bonds and potentials are imposed on a uniform 
hopping model with t = 1.0, and for solid line, the amplitudes 
of disorders are St — 0.2 and Sv — 0.2; while for the dash line, 
the amplitudes of disorders are St = 0.3 and Sv = 0.3. One 
sees that slight disorders can induce a peak-like structure at 
zero energy. 

hard-wall boundary point. In this case, for any set of 
{ti}, the non-Hermitian matrix A that enters Eq.@ is 
given by 

/ti -\ 

For the zero-energy state, because A</> e = 0, it is easy 
to see that 4> e — 0; while the wavefunction at odd 
sites is determined by A tp = 0. If we set \I/q = 1 , 
the wavefunction at site 2N + 1 is given by ^q N+1 sa 
(t 1 t 3 t 5 ...t 2 N-i)/(t 2 t4:t 6 ...t2N)- Let Xi = \t 2 i-i/t 2 i\, then 
lnl*^ 1 ! sa Etihi^- Clearly, because ln|*^ +1 | - 
In l^o^ ^l ~ hi x n , the logarithm of the wavefunction at 
odd sites behaves effectively as a random walker. Since 
a random walker has large probability to go to ±00, \I/o 
has high probability of decaying to zero far from the edge 
and becomes a localized zero-energy mode. This analogy 
leads to ((In \^l N+1 \) 2 ) = VN<J with a being the stan- 
dard deviation of Inx. In other words, l*^ 1 ! ~ e ±a ^ 
where ± correspond to states localized at either ends, 
respectively. We emphasize that this analysis indicates 
that only the standard deviation of In x is relevant and 
there is no need for the assumption of Gaussian-type 
randomness, which is often invoked in previous works. 
Furthermore, the random-walk nature makes the zero- 
energy peak much more easily formed. This is illustrated 
in Fig. 0] where we show that even slight disorders near 
the edge may induce features resembling zero-bias peaks 
in tunneling measurements. In this case, the zero-energy 
state tends to decay from the edge but will not become 
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(a) 




where 




(b) 




FIG. 5: Schematic plots of (a) the spectrum with a midgap 
state and (b) the spectrum with an impurity state, which 
may appear at the position of the solid line or the dashed 
line. Clearly, (b) may be viewed as the square of (a) (see text 
for details), i.e., if one folds (a) with respect to the midgap 
state, one gets (b). 



localized. Instead, after joining the non-disorder bulk 
region, it becomes a resonant stated. Such phenomena 
may have already been seen in experiments— 

IV. IMPURITY STATES AS MIDGAP STATES 

In this section, we demonstrate that in the supersym- 
metric approach the ordinarily recognized impurity states 
can be viewed as disguised midgap states. Let us con- 
sider the Goodwin model for the surface state—. In this 
model, it was proposed that the surface state arises be- 
cause the potential suddenly changes near the surface or 
the edge. In the tight binding limit, the Hamiltonian is 
given bjiSi 



i=l 



tc]cj_|_i + h.c. + Uc\ci. 



(23) 



In other words, there is an impurity potential localized 
at the first site. It is commonly recognized that under 
appropriate condition, an impurity state (in this case, it 
is the Goodwin edge state.) may arise and exhibit as an 
isolated line in the spectrum as illustrated in Fig. Efb). 
Clearly, we see from Fig. \5\ the spectrum with an impu- 
rity state is essentially the square of spectrum shown in 
Fig. Ela), i.e., the spectrum with a midgap state. 

To establish the relation described above, we go back 
to previous analysis on the ti-t2 model. As is obvious, 
the spectrum of H s (= (H2) 2 ) is the square of that for 
H2, fulfilling the relation indicated in Fig. [SJ It is hence 
useful and more transparent if we explicitly write down 



AA* 






A f A 



(24) 



/ A 



iff = AA f = 



ht 2 



\ 



t\t2 t\ + t 2 



q + t 2 



(25) 



is the effective Hamiltonian for the odd sites and 



iff = A f A = 



tf+ti 




ht 2 
t? + ti 





t 2 + t 2 



(26) 



is the effective Hamiltonian for the even sites. One sees 
that iff and iff differ by the potential at site 1 . As 
mentioned, this is entirely due to the missing bond cut 
off by the boundary. On the other hand, even though 
H s is block-diagonal, it does not imply even and odd 
sites are independently from each other. In fact, they 
are connected by the supercharge H 2 - The point is that 
except for the zero energy which is an eigenvalue of iff , 
iff and iff share the same energy eigenvalue E 2 with 
E being the spectrum of if2. The supersymmetric rela- 
tion between even and odd sites enables one to solve the 
Goodwin model Hq as follows. One first rewrites 



E 



[(ht 2 )<H 



+ h.c. + (tl+1%)c i Ci 



t\c\c\. 
(27) 

Clearly, it shows that the effective Hamiltonian for the 
odd sites ie equivalent to the Goodwin model with t = 
t 1 t 2 , U = —t 2 , and = — (if + t 2 ). The wavefunction of 
the midgap state on odd sites then become the wavefunc- 
tion of the impurity state in the Goodwin model. Fur- 
thermore, the energy of the impurity state can be easily 
found to be Ei m = — [t\ + if). By solving t\ and t 2 in 
terms of t and U, we find Ei m = U + t 2 /U. Since the 
midgap states exist when < ^l, the impurity state 
appears only when \t\ < \U\, in consistent with the stan- 
dard approach— The discussion above concerns with the 
case U < 0, hence the impurity energy resides on the left 
side (the solid line) in Fig.|Sfb). For U > 0, the Goodwin 
model maps to —iff. One obtains the same expression 
for the impurity energy Ei m — U + t 2 /U except that it 
resides on the right side (the dashed line) in Fig. [5£b). 

In addition to the energy of the impurity states, the 
above analysis also implies that the entire spectrum is 
simply Eq = 2i cos ka. Furthermore, the supersymmet- 
ric relation between even and odd sites enables one to 
write down all wavefunctions for the Goodwin model ex- 
plicitly. This is entirely due to the fact that the Hamil- 
tonian of the supersymmetric partner to the Goodwin 
model is iff which is a uniform hopping model. We ob- 
tain the wavefunction for the impurity state ^ci 71 ) = 
(t/U) n ^ 1 , while for the extended states, when U < 0, 
iF G (n) = [y/\U\ s'mnka + t/y/\U\ sin(n - l)ka]/E k with 
E k = ±^\U\ + t 2 /\U\ + 2t cos ka, and when U > 0, 
* G (n) = [\f\U~\ sinnka - t/y/\D) sm(n - l)ka]/E k with 
E k = ±y/\U\ +t 2 /\U\ - 2£ cos kci. 
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(a) Soliton 




(b) Impurity 



t'l-t^-t's-t'^ model, i.e., the p — 4 (H4) model. Specifically, 
one finds that Hq is equivalent to the block in H 2 that 
describes the odd sites; while its supersymmetric partner 
is the original t\-t% model with no impurity potential. 
The impurity state of Hq is identified as the midgap state 
of the H4 model at E — 0. Simple calculation then yields 
E lmp = h 2 /(2U) + U/2 + ^/(t 2 2 /(2U) + U/2) 2 + t 1 2 . 
Thus, one concludes that in general the impurity state 
of H p can be mapped to the midgap state of H2 P - 



u 



FIG. 6: The equivalence of the electronic state in a SSH 
soliton background (a) and in an impurity background (b). 
The electronic state for an impurity resides only on the black 
points. In (a), the thin line represents hopping amplitude ti, 
while the thick line represents hopping amplitude £2- In (b), 
the hopping amplitude is t\t2. 



One can also apply the same approach to the bulk 
case. In Fig. Efa), we show that the soliton con- 
figuration in polyacetylene discussed by Su, Schrieffer, 
and Heeger (the SSH soliton)iiii is essentially a bulk 
impurity illustrated in Fig. Elb). The same analysis 
leads to the results: 



and * lmp 
for U > 



(n) 
0, 

\U\/(2t) + ^(UjW 



For U < 0, E imp 
\U\/(2t) — 



E, 



imp 



= VTP 

M 



VW/W+ 

+ At 2 , and 



VIP 

\n\ 



At 2 
while 



imp 



(») 



1 



. Here n is measured from 



the impurity site. Note that unlike the case for semi- 
infinite chains, the impurity state for the bulk ID chain 
exists for any U and t. 

The purpose of the above analysis is only for demon- 
stration. In principle, it can be also applied to cases when 
the potential extends from the first to other lattice points 
in the Goodwin model. For instance, one can include an 
additional term Vc\c2 in the Goodwin model. In that 
case, the first and second bonds in the corresponding t\- 
t 2 model has to be different and denoted as t\ and t' 2 . One 



t\ 2 



ti 



to 



then obtains that t — t\t 2 = t'it'2, U 
and V = t' 2 — t2 2 ■ Since for any changes of finite num- 
ber of hoppings, the energy of the edge mode stays at 
zero. The energy of the impurity state is still given by 
Ei m p — — (ti 2 + i 2 2 )- Solving t\ 2 and t 2 2 in terms of t, 
U, and V yields the energy of the impurity state. 

We close this section by pointing out that the idea of 
mapping impurity states to midgap states is quite gen- 
eral. In addition to the simplest version of the Goodwin 
model discussed so far, it also works for more general 
cases. For instance, consider the generalized Goodwin 
model Hq in which one introduces an impurity poten- 
tial at the 1st site of the t\-t2 model. This model can be 
mapped to the effective Hamiltonian for the odd sites of a 



V. SUMMARY AND OUTLOOK 

In summary, in this work we have shown that the prop- 
erties of midgap states in semi-infinite ID nanowires are 
dictated by a underlying discrete supersymmetry. This 
supersymmetric structure generalizes the ordinary super- 
symmetric quantum mechanics and offers a new point of 
view toward the origin of edge states. In the presence 
of hard- wall boundary condition (-0 = 0), the sublat- 
tice which directly connects to the hard-wall spans the 
null space of the supersymmetric ground state. As a 
consequence, the energies of the midgap states are de- 
termined by the eigenvalues of a reduced Hamiltonian, 
Ea. (|2L)|) . whose size is much smaller than that of the 
original Hamiltonian. This reduction in matrix size sig- 
nificantly reduces the computation cost for determining 
the occurrence of the midgap states and offers a way to 
manipulate them. As an application, we investigate a 
structure with superlattice in hopping. In this case, the 
number of edge states is simply the period of the super- 
lattice minus one. Therefore, changing the period offers 
a way to control the number of localized electronic states 
at the edge of the nanowires. 

While so far in this work we have not considered 
electron-electron interactions, from adiabatic continuity, 
the results obtained here should still hold for cases when 
the electron-electron interactions are weak (so that quasi- 
particles are well defined). Moreover, if the "ID chain" 
results from the reduction of a higher dimensional struc- 
ture where interactions are not important, then it is le- 
gitimate to ignore interactions in its effective ID model. 
This is of course not correct in truly ID atomic chains 
where it is known that interactions may dominate the 
physics and the quasi-particle pictures may fail. In this 
case, however, the states we obtain can be used as the ba- 
sis to express the full Hamiltonian (with interactions) uti- 
lizing the relation c ia = i/j (i)c 0a + Y J E^E{i)c E(T . When 
Coulomb interaction is included, it reduces to the An- 
derson model 2 ^ in which the edge state acts as an im- 
purity state. The scattering of extended states by the 
edge states essentially causes the Kondo effect, result- 
ing in the zero-bias peaks near the Fermi energies*. On 
the other hand, the interaction also correct the localized 
edge state. This is conventionally analyzed in the Fano- 
Anderson model 22 , in which it is known that as long as 
the new energy found remains inside the gap, the corre- 
sponding state is a localized state. In either of the above 
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mentioned effects, our results will serve as useful inputs 
for attaining the final corrections. 



Ea. (|18|) . It is then suffice to prove [Hs, q P ] 
that Hs has the following generic form 



0. We note 
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APPENDIX A: SUPERALGEBRA IN INFINITE 
p-PARTITE SYSTEMS 



Ha 



/Sn 
S21 
S31 



S12 
S22 
S32 



S13 
S23 
S33 



\ 



-i,p- 




1 





0/ 



-H° s 



H 1 



Here H s 



(Bl) 

has the same form as that of the SUSY Hamil- 
tonian for the corresponding infinite chain except that 
semi-infinite lattice points are removed, hence it is block- 
diagonal with the form: (Hg) nm = So5 nm , where n and 
to are the block indices and 



In this appendix, detailed superalgebra behind our 
generalized SUSY quatum mechanics is presented. For 
an infinite p-partite system, after modulo p, the system 
reduces to the set {1, 2, 3, • • • ,p\ with each number rep- 
resenting different sublattices. The reduced system is 
periodic with p + 1 being identified as 1. In this periodic 
space, a set of generators {Q n m', n,m = l, 2, ...,p} can be 
defined. Here Qnm are p x p Hermitian matrices whose 
only nonvanishing elements are 1 in the n th row and m 
th column and the m th row and n th column. Note that 
Qnn has only one element, 1, in the nth element along 
the diagonal. Obviously, when n ^ to, Q nm permutes 
the subwavefuncions </> n and (j> m ; when combined with 
the hopping strength t nm , in addition to permutation, it 
also rescales the wavef unctions. The Lie algebra formed 
by Qnm is a superalgebra because the anticommutator 
is necessary in order to be closed. The followings are 
nontrivial commutation relations: {Qi m ,Qmn} = Qin 

for I ^ Tl , {Q nm , Qnm} = 2Qnn + 2Qmm for 71 ^ TO, 

and [QnmQmm] = SnmQn', ah the other commutators 
are zero. It is straightforward to check that for infinite 
systems, the SUSY Hamiltonian Hs defined in Sec. Ill 
commutes with all Q nm even if the operation of modulo 
p is not performed. 



APPENDIX B: DERIVATION OF [H S ,Q V 



In this appendix, we outline the proof of [Hs, Q p ] = 
for semi-infinite chains. We first write Q p = q p +q P where 
q p is obtained by setting the last column in Q p to zero in 



/ t •• \ 

t t ■■ 

t •• • (B2) 
t 

\ o o •• t o / 



and t — £^2^3 ■ -t p . The block matrix S nm represents 
missing hopping amplitudes that goes between sublat- 
tices m and n due to the presence of the boundary point 
at i = 0. When computing [Hs,q p ], one needs to com- 



pute [S , A lp ], A lp • Si m and A p _ l p • S p _i, m , thus only 

Si m and S p _i. m are needed. It is straightforward to show 
t 

that [So, A lp ] has only one element 
element along the diagonal (= 5±x) 



-tt p , which is the 1st 
To obtain S nm , one 



needs to multiply H p to itself by n (< p) times because 
Hs is a polynomial of H p containing at most pth power 
of H p . Now the multiplication of H p to itself n times ef- 
fectively hops a particle n times. Since n < p, only when 
the particle starts from the lattice points 1 < i < p — 1, 
it will have chance to visit i = and thus will have miss- 
ing paths when i — is removed. It implies that all S nm 
have only one element, which is also the 1st element along 
the diagonal. In addition, S p _i jITl = for to > 2. As a 

result, by using Ea. ((TC|l . we obtain A lp • S lm = and 

hence we only need to compute Sp-i^i. Since the only 

missing path between the point 1 and the point p — 1 

isl— »0— ► 1 — ► 2 ■ • • — * p — 1, we find that S p _i ; i 

is ~t p 2 tit2 ■ ■ -t p _i(5ii. This result, when combined with 
t 

[Sq, A lp ] = —ttpSn, we finally obtain [Hs, Q P ] = 0. 
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